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Abstract 

We consider a two-dimensional singularly perturbed transmission problem with two 
different diffusion coefficients, in a domain with smooth (analytic) boundary. The so- 
lution will contain boundary layers only in the part of the domain where the diffusion 
coefficient is high and interface layers along the interface. Utilizing existing and newly 
derived regularity results for the exact solution, we design a robust hp finite element 
method for its approximation. Under the assumption of analytic input data, we show 
that the method converges at an exponential rate, provided the mesh and polynomial de- 
gree distribution are chosen appropriately. Numerical results illustrating our theoretical 
findings are also included. 
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1 Introduction 



The approximation of singularly perturbed problems has retained the attention of many 
authors in recent years. Let us mention [51 El EDI Ell E31 El] and the references quoted there. 
However, in all references quoted no analysis is carried out for differential operators with 
piecewise constant or piecewise smooth coefficients. On the other hand, in many real life 
applications, the differential operators have such piecewise coefficients that may have a very 
large discrepancy. In that case, the solution of the problem will contain boundary layers near 
the exterior boundary (as usual) but will also contain interface layers along the interface where 
the coefficients have a large jump. We refer to [I] for the description of this phenomenon in one 
and two dimensions and to [12] for several numerical methods for the robust approximation 
of such problems in one-dimension. 

The goal of the present paper is to extend certain results from [12] to two- dimensions. In 
particular, we consider a singularly perturbed transmission problem in a domain with analytic 
boundary. Under the assumption of the data also being analytic, we provide an asymptotic 
expansion for the solution (in the style of [6]) that provides the necessary information for the 
design of a robust finite element method that converges at an exponential rate as the degree 
p of the approximating polynomials is increased. The expansion of the solution includes 
an outer (smooth) part, an inner (boundary layer) part, an interface layer and a (smooth) 
remainder. The regularity of each compoment is studied and known results from [8] allow 
us to treat the outer and inner parts, as well as the remainder (defined on one part of the 
domain). The results obtained for the regularity of the interface layer (and the remainder 
defined on the other part of the domain) are new and in line with those reported in [12] for 
the one-dimensional analog of our model problem. Our work closely follows what was done 
in [7] but also includes the additional analysis for the interface layer. 

The paper is organized as follows: In Section [2] we present the singularly perturbed problem 
and describe the typical phenomena. Section [3] is devoted to the expansion of the solution of 
our model problem into the parts mentioned above (i.e. outer, inner, interface and remainder). 
The regularity of each component is also described in that section. Section 2] gives the 
main approximation result and in Section [5] we show the results of numerical computations 
illustrating our theoretical findings. We end with some conclusions in Section [6j 

Throughout the paper the spaces H S (Q), with s > 0, are the standard Sobolev spaces on 
the domain Q C M 2 , with norm || • || Sj q and semi-norm | • | S) n. The space Hq(Q) is defined, 
as usual, by Hq(Q) := {v G H 1 ^) : v\ dn = 0}. L P (Q), p > 1, are the usual Lebesgue spaces 
with norm || • ||o, p ,n ( we drop the index p for p = 2). Finally, the notation A < B means 
the existence of a positive constant C, which is independent of the quantities A and B under 
consideration and of the parameter e, such that A < CB. 
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2 The model problem 



Let f2 + and f2_ be smooth domains in M 2 , with respective boundaries dQ + and <9f2_, such 
that dQ + fl <9f2_ = S; an example is shown in Figure 1 below. We assume that dQ is an 
analytic curve, i.e. dQ± and £ are analytic curves. Moreover, we assume that dQ + \E, as well 
as £ are connected. We will write f2 = Q + U and for any function u defined on Vt we will 
denote by u + (resp. uJ) the restriction of u to f2 + (resp. VtJ) and we will write u = (u + , uJ). 




Figure 1: Example of the domains Q + and 
We consider the following singularly perturbed transmission problem: Find u £ = (n 



£ 1 £ 



such that 
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where A denotes the Laplacian operator, e G (0,1] is a given parameter, f±,h are given 
smooth functions and v denotes the outward normal vector along S oriented outside Q + . The 
formal limit problem of ([I))-©, as e — > 0, is 

u° + = f+ in n+, 

-Au°_+u°_ = /_ in Q_, 

u° + = on 

u°_ = 0on9n„\E, 

u° + -u°_ = on E, 
du° 



Since, in general, /+ does not satisfy the boundary conditions / + = vP + on dfl + \L and 
/+ = u°_ on S, we expect that the solution u e will contain boundary layers along dQ + \E and 
an interface layer along S. 

We assume that the data of our problem is analytic and satisfies 

(7) \\^ P f±\L,n ± < C/ ± T? ± P! V p = 0, 1, 2, 

(8) l|V^IL iE <C fc7 £p!Vp = 0,1,2,..., 

for some positive constants C/ ± , 7j ± , C^, 7 fe , where Vs denotes the tangential derivative along 
S. The following theorem gives bounds on the derivatives of the solution to (JTJ) — (JSJ) that are 
explicit in terms of the order of differentiation as well as the singular perturbation parameter 
e. 

Theorem 1 Let u e = be the solution to (CZP-(0|) with the data satysfying (0), (Q|). 

Then there are constants C, K > depending only of the data such that 

(9) £||D a ^|| on + ||^ a ui|| on _ < CeK |Q| max{|a| Va = 1,2,.... 

Proof. This follows from the local estimates 

£\u+\2,B X0 nn + + \u e _\2,B X0 r\Si- < C(^ 1 \\f+\\o,B' xo nn+ 
+ ||/-||o 1 fl4 nn_ + ||/i||Bi ns3 + e||w+||i,Bi nn + + ||«I||i,B4 nn_), 

for all sufficiently small balls B xo C B 1 centred at Xq G S (proved by a local change of 
variables and some reflexions to reduce the transmission problem into a Dirichlet problem 
and a Neumann one in half-balls) and the use of Morrey-Nirenberg techniques (see Theorem 
2.1 in [7] or Theorems 5.2.2 and 5.3.8 in [3]). ■ 

It should be noted that (Q gives sufficient information for the approximation to u e in the 
so-called asymptotic case, i.e. when the degree p of the approximating polynomials satisfies 
p > 0{£^ v ). For the pre-asymptotic case, i.e. when p < 0(e _1 ), we will need the regularity 
results provided in the next section. 

3 Expansion of the solution 

The solution of ([I])-© may be decomposed as 

(10) u £ = w £ + XblU BL + Xi L u £ IL + r\ 
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where w £ denotes the outer (smooth) part, u e BL denotes the boundary layer along <9f2 + \E, 
u s IL denotes the interface layer along E and r e denotes the remainder. The functions XbliXil 
denote smooth cut-off functions (see equations (ITTj) . ffTgj) ahead) in order to account for the 
fact that the aforementioned components do not have support in the entire domain Q. 

In order to define the inner (boundary layer) expansion we introduce boundary fitted 
coordinates as follows: Let (X(6), Y(9)) ,9 G [0, L] be an analytic L— periodic parametrization 
of <9f2 + \E (by arc length), such that the normal vector (—Y'(9),X'(9)) always points into the 
domain Q + . Let k+(9) denote the curvature of <9f2 + \E and denote by Tx, the one-dimensional 
torus of length L. By the analyticity of dQ we have that the functions X, Y and k are analytic. 
We also let p > be a fixed constant satisfying 

(11) 0<p < 



\ K +\\L°°({0,L)) 



Then the mapping ip : [0, p ] x T L — > Q + given by 

(12) J, : (p, 9) (X(9) - pY'(9), Y{9) + P X'{9)) 

is real analytic on [0, p ] x Xl- The function ip maps the rectangle (0,p ) x [0, L) onto a 
half-tubular neighborhood VL° + of dQ + \E, which may be described as 

(13) Q° + = {z - pn z : z G <9fi+\E, < p < p } , 

with z = z{9) = (X(9), Y{9)) and n z the outward unit normal at z G dfl + \T,. 

The interface layer will also be defined in a neighborhood of the interface E. Quite 
analogously, let (X^(9),Ye(9)) ,9 G [0, Ls] be an analytic L^— periodic parametrization of E 
(as above), let kt,(9) denote the curvature of E and denote by Tl e the one- dimensional torus 
of length Ls- With p E > a fixed constant satisfying 

(14) < Ps < 



iI k £|| LO o ([0iLe)) 
we define, analogously to ( TT31 . 

(15) til = {2-pn E : z G E,0 < p < p E } , 

with z = z(9) = (Xj^(9), Yz(9)) and n E the outward unit normal at z G E. 

The smooth cut-off functions x BL , Xil appearing in (|T0|) are defined as follows: Let p x , p 2 
be given satisfying 

(16) < Pl < Po , < p 2 < p E 
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and let Xbli Xil be denned on f2 + via 



(17) XbM 



1 for < dist(z, <9fi+\E) < Pl 

for dist(x, dQ + \E) > {p x + p ) /2 



1 for < dist(x, E) < p 2 
for dist(x,£) > (p 2 + p s )/2 



(18) Xii(z) = 

The above will be utilized in sections 13.21 and 13.31 ahead. 



3.1 Construction and regularity of the outer part 

We begin by constructing the outer part w e in ( ITOj) . To this end, we expand the solution 
u e = (u^_, u £ _) as a formal series in powers of e, 

(19) u% = + euf + e 2 uf + ... 

and insert it in the differential equations (JTJ-flSJ), equating like powers of e. This allows us to 
get expressions for the functions uf,j = 0, 1, 2, .... In particular, for ul- we obtain 

(20) u+ = /+, u+ = A^' +2 )/ + , = 0, j = 1, 2, ... 
where A^ denotes the iterated Laplacian. For Uq we obtain 

(21) - Au + = /_ in 

(22) «o = 0on<9fi_\E, 

(23) ^ = -ZionE. 

ov 

For j > 1 we find w^,_! = and 

(24) -A«-+«2j = 0infi_, 

(25) u 2j = 0on9n_\S, 

(26) &7 = ^7- onS - 

Note that u 2 j is not explicitly known but is solution of a Dirichlet-Neumann problem in 
Due to the analyticity assumption, u 2 - is analytic as well (see equation ( 1341) ahead). 

Using the above, we can define the outer expansion as 



M 

(27) w ± = w± = J^fc 

j=0 



where M is the order of the expansion (i.e. the number of terms that we will include) and 
will ultimately be taken to be proportional to 1/e (cf. [5], [7]). It is not difficult to see that 



(28) 
and 
(29) 



(-e 2 Aw 



M T % 



) - /+ = 5 2M+2 A (A/ + 1) ^ 



("A 



w 



A I 



W 



AI 



)"/- 



Moreover, we have the following theorem. 



Theorem 2 Let w M be defined by pT7| ). T/ien i/iere ernst positive constants K\ and C de- 
pending only on the data of the problem, such that if eM is sufficiently small then 



(30) 



\]J a w + II <" &H a ll~.li W „. r- rj2 



M|| o ! Q + 



< Kp\a\\Vae N 



o- 



(31) 



Proof. From Theorem 2.2 of [8] we have that 

ll^<IL,n + £ A f (l + (2MeK 2 ) 2M ) V a e N%, 

so if 2MeK2 < 1 we get (1301) . In order to establish ()31| . we first consider u$, which satisfies 
the Dirichlet-Neumann problem (l2~TT)-(l23l). Since the data of this problem are analytic, we 
have that Uq is also analytic [3J, and moreover 



a 



o \k,U. 



< C k+l k\ V k e N . 



Next, we consider u 2 j, j = 0, 1, defined by (j2"lj) - (j26|) . with again the data being analytic. 
Casting (|24j) - (j26p into a variational formulation, allows us to write 



it 



2 illi,n_ 



dutj-2 - / 



<9z/ 



r2jlli,n_ ' 



0,2 



which, using ([20]) . gives 
(32) 



M 2j|lin ~ 



dv 



< 



0,£ 



a(A 2 V-t 



<9z/ 



< l|v 4j+ V- 



o.s 



From [3] we have that there exists C G 1R + such that 



(33) 



fe+i 



9m 



2j 



and we note that (see eq. ( |26l) 

2 



du 



2j 



du 



du 



2j-2 



du 



du 



< ||A 2 V4 



2 

\e+i,n + 



< 



\a\<t+l ^ 

Hence, (1331) becomes (with the aid of ([7]) and ([32 

A- 



£ J\D a A^f + \ 2 dx< \\D a A 23 f- 



2 

^~ 1 1 oo,f2j 



|«|<£+1 



fc\ I 2 ilfc,Q_ 



* Et[ E II^ 2 v+IU + + I|a 2 V- 



<£\ 

=1 |o|<^+l 



+ lloo,E 



S Cfc+1 'E^[ E 7 H+2 MI«|+2j)! + 7 2 M2j)! 

=1 '|a|<f+l 



< 



C fc+1 ^V +2j (^ + 2j)! 



< C fe+1 fc 2 (2j)!^7 fc+2j f 



2j 



< C 1 fe+1 (2j)!(l + 7) fc+2J 7 2j 

< C 1 fe+1 (2j)!7i J 

for a suitable Ci, 7 X > 0. This shows that u^a are analytic and V j = 0, 1, ... 



(34) 



u 2j\ kn <C k 1 +l k\{2 ] )\^ , ken. 



Thus, from the definition of w M we have 



M 



\W 



2.; 



j=0 j=0 



< 



,2/ 



Cj+^JV* (2Mf 7 2j < C^kl^, (2^M 7i ; 
i=o i=o 

< C k+1 k\, 

provided 2eM 7l < 1 (so that the above sum can be estimated by a converging geometric 
series). Estimate f[3"Tj) follows. ■ 

Remark 1 The above theorem gives bounds on the smooth (outer) part of the solution to 
(QP-(0|) under the assumption that eM is sufficiently small. In the complementary case, the 
asymptotic expansion loses its meaning. 
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3.2 Construction and regularity of the boundary layers along dQ + \E 

Boundary layers are introduced in order to account for the fact that the function does not 
satisfy the boundary condition on <9f2 + \E (cf. (1281)). These are precisely the ones constructed 
and analyzed in [8], so we will only outline the procedure and quote the relevant results from 
[8]. The boundary layer correction u e BL of is defined as the solution of 

(35) L e u £ BL = OmQ + , 

(36) u% L = -wjj on <9fi+\£, 

where L £ is defined as 

(37) L e u := —e 2 Au + u. 
With k+(6) the curvature of <9f2 + \£ we set 



and we have (see, e.g. pQ) 

Au(p, 9) = d 2 p u - n + (9)a + (p, 6)d p u + a\(p, 6)d 2 e u + P k' + (6)(j%{ P , 9)d e u. 
Introducing the stretched variable p = p/e, the operator L £ becomes 

(38) L £ = -of + Id + eK + (9)a + (ep, 9)ch p - e 2 a\{ep, 9) - e 3 pK + (9)a 3 + (ep, 9)d e . 
Expanding the above in power series of e, we can formally write 

oo 

(39) L £ = 

i=0 

where the operators Lj have the form (see equations (2.12)-(2.14) in [8]) 

(40) L = -of + Id, Li = -p^af'dp ~ P^'d 2 - p l ~ 2 al 3 d e , i > 1, 
and the coefficients a* are given by 

(41) a\ = - [n + (9)] l+1 , 4 = (* + !) [*+(*)]' , 4 = { * + ^ + 2) [«+(*)]< i G N , 



(42) a\ = 4 = 4 = for i < 0. 
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We next make the formal ansatz 

oo 

and insert it into f )35|) . This yields 



(43) Y/T, L Pi-i = °» 

i=0 j=0 

allowing us to find the following problem for the functions Ui(p,9),i = 0, 1, 2, ...: 

(44) - d$Ui + Ui = Fi =: Fl + Ff + Ff, 

i— 1 i— 2 i— 3 

(45) I? = Y?*idfii-i-k » *? = 5>*^ft-2-* , i? = ^P fe+1 a^ e ^_ 3 _ fc , 

fc=0 fc=0 fe=0 

where empty sums are assumed to be zero. (See, also, equations (2.15)-(2.16) in [8]). The 
above are supplemented with boundary conditions 

Ui — > as p — 7- oo, 

-[A^ 2 )/] anxs if zGNiseven, 
if z G N is odd. 



dfl\Y: 



The boundary layer (inner) expansion in (TTUT) is then defined as 

2A/+1 2M+1 

(46) u £ BL = < L (p, 9) = £j Uj(P, 0) = E £ ^G%, 

j'=o i=o 

and by construction, it satisfies the boundary condition 

2A/+1 

[ u sJan\E = ~ E ^ [ AW ^an\s ■ 

i=0 

By Theorem 2.2 of [8] we have that for every a G [0, 1) and all p, m G Mo, 

2M+1 N 

(47) 



for 9 G Tl, p G [0, p ], with i^i, > independent of e,p and m. Moreover, by Lemma 2.12 
of [Sj, there exist constants K, > independent of e such that 

(48) \L £ u A J L (p, 9) I < K 2M+2 (e(2M + 2) + |p|) 2M+2 V (p, 0) G B Po (0) x 5(6), 
where -85(2;) denotes the (open) disc in the complex plane of radius 5 centered at z, and 

(49) S{Q) = {9eC: lm(0) < 6} . 
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3.3 Construction and regularity of the interface layer on E 

For a function oj = (u + ,U-) we denote the jump [[u;]] s on E as 

(50) [H] E := K)| s - 

We define the function t>/(p, #) := (vj, vf) as the solution of the following problem: 

-e 2 Avf + vf = in fi + 
— Ai>J + = in Q_ 

oo 

(51) { N^-E^Ws 

j=0 

2 ^ _ <9"f \ 

With p = p/e as before, we write 



i=o 



9«+_ _ au^ 

9p 9p 



S+(p,0)=<;+(p,0), 

and problem (15T1) becomes 

f (-of + Id) + e« + (0)<7(ep, 0)6>£+ - e 2 a 2 (ep, 9)vf- 
-e 3 pK' + (9)a 3 (ep,6)devf = in fi+ 
— Ai/^ + f ^ = in Q_ 

oo 

( 52 ) S - vj) | s = [K-]] s 

i=0 



9p 9p 



j=0 



9p dp 



Now, we write 
(53) 



v 



t = , vj = ^ j Vf, 

3=0 j=0 

and insert it in ( |52|) equating like powers of e, to get (utilizing again the expansion (|39|) ^ 

dft> + + i? + = F/ + Ff + Ff in R + V j > 
AVf + Vf = in fi_ V j > 



(54) 



(uj - U;y) on E V j > 



^t+i-^+i) = 0onEVj>0 



At/ - — - a -T/+ 

5p K 2j &p V 2j-l 



aV ~ _ du~ 
dp dp 

9 -u~ 



on E 

9 „,+ 



9p "2j a p <*2j- 



<-_ 2 onEVj'M 



At/- _ At7+ 

9p V 2j+1 dp V 2j 



on E V j > 
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with Fh F?, Fj given by ( 145|) but with U replaced by V + . So for j = 0, we have 



(55) { "t = ?onE 



-AV ~ + V ~ = in fi. 

av o~ _ duo 
dp dp 

Vq = on <9Q_\E 



(56) 



-of F + + V = in K+ 

K) + = - (*4 - u o ) 011 s 



-AVf + 1/f = in a 



(57) { %K = on E 

Vf = on <9fi_\E 



(58) 

In general, for j > odd we have 



-of Vi + + t>+ = t? + in 
7+ = on E 



-AV 2 ~ +1 + V 2j+1 = in n. 



(59) < |^ +1 = |Kj on S , 

= on <9fi_\E 



(60) 

and for j > even we have 



V 2j+1 = V 2j+1 011 S 



-AV£ + V 2j = in n. 



< ( 1^ - = - (|s- - 14- 2 j on s > 

V 2j = on <9fi_\E 



(62) 



-<9jVg + V = F 2 ) + F|. + in 

V 2j = V 2j ~ ( U 2j ~ U 2j) 011 S 



The regularity of the functions Vj , is given by Theorem H] below. For its proof, we 
will need the following lemma. 
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Lemma 3 Let Uj (p, 9) ,j = 0, 1, 2, be the solutions to 
(63) 



-d$Uj + Uj = Fj (p, 6) in 
Uj = Gj (9) on S 



where Fj (p, 9) = Fj (p, 9) + (p, 0) + F? (p, 0) owen fry (g^-(g3]] and Gj satisfy 

(64) |Gj(0)|<Cb7^, 

for some positive constants Cg,1g depending only on the data. Then, there exist positive 
constants Q^Cu^u, depending only on the data, such that 

(65) \Uj (p, 9) | < Cut*, (1 + j + pf e"? V (p, «)Gl + x5 (6) , 

where S(Q) is given by (fT^P - Moreover, for any a G [0, 1) i/iere exists K G M + depending 
only on the data, such that 



(66) 
and 



|^(p^)|<A^y(l-«)" J e-^, 



(67) \d^dpj (p/e, 9) | < ^V 1 "^ (p + 1) 1/2 g! (2/0f 7^ (1 - e*p e"^ V p, g G N„. 

Proof. This is essentially a combination of Lemmas 2.9 and 2.11 in jH]. Estimate (1651) 
follows directly from Lemma 2.11 in [8], while fl66|) follows from ( 165]) and Lemma 2.8 in [8]. 
Finally, ( 16T|) follows from Cauchy's integral formula, in exactly the same way as in the proof 
of (2.24) in [8]. ■ 



Theorem 4 Let V, satisfy (C2|) ; / fgi]) and satisfy [6tyl , / fg^j) . Then there exist constants 
C, 7, K, O > depending only on the data, such that 



(68) 



u>/w/e /or 6* G Tl e , p G [0, p E ], a G [0, T 



(69) 
and 
(70) 



K + (p/M) <^f(l-a)- J e 



<^V; + (p/ £ , f) < e~ p e^ (p + 1) 1/2 q\ (2/0)* i^f (1 



forp,q G No and # G 5(6) ou>en 6y ( f^P - 
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Proof. The proof is by induction on j. First we note that estimates ( |69|) . (!70|) follow from 
Lemma [31 provided we show that flM|) is satisfied, i.e. on E the functions are bounded by 
Cji^i for suitable constants C, 7 > 0. This will be verified during our induction argument; in 
fact it will be the only thing we will show for Vj + , with the understanding that an application 
of Lemma [3] gives the desired result. 



For j = we see from the variational formulation of ( )55|) that ||V || < 
hence by ( |34"|) and [3] 

1 



u 



Mi o_ ' 



fc+i 



fc+1 




II 








Next, for V + we see from fl56l) that Vq + (p, 0) = Gq{Q) e~ p l £ for some function Go (9) that 
depends on V^~,Uo,«o- By the above, (120)) and (1341) . we have that in the case j = 0, the 
boundary data for is bounded by Cjirf for suitable constants C, 7 > 0, hence by Lemma 
[3] the bounds flTD]) hold for V+. 

Now, from the variational formulation of (l57j) and the fact that V^ + (p, #) = Go(9)e~~ p , we 
have |^ + (p,#) = -G o (0)e^ and then ^V + (0,9) = -G (9), hence 



where 
(71) 

Thus, 



/ ( Wf ■ W + V^~V) dx = - [ G V dx V F e F* 1 ^- 
= { M G i/ 1 (n_) : ula^j, = 0} . 



< 


l|C 


o,s ~~ 


V + 


0,E 


K 


< 




lo,E + 


\ u o\ 


0,S 


u o\\ 
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From [3] and the above result, we get 



k\ 



lifc,Q_ 



< c 



k+l J 



k-2 ^ 
fc-2 



dv n - 



< c k+1 

< c k+1 



dp 

£=0 
fc-2 1 



Ci 



k 



+ » 



N<+i,E 



+ ||«o IL+|,s) + ^ 



< c 



fc+1 



which leads to 
(72) 



\k,Q- 



< C k+1 k\. 



In an analogous way as V^(p, 6) = G (6) e~ p , we find that Vf(p, 9) = -$a%(9)e-P+Gi (6) e~P, 
for some function G\ (8) that depends on V{~ , hence in view of (172)) we see that the boundary 
data for satisfy the appropriate bound. As a result, (I55|l - (l70l) hold for V7 1- as well. 

So, we assume that fl68|) - (j70l) hold for j and we will establish them for j + 1. 

The case of odd j: If j is odd, then j + 1 is even and we would like to establish bounds 
for and (with 2s = j + 1), which satisfy (I6TI) . (162|) respectively. First, for V 2 ~ we see 
from the variational formulation of (l6Tj) that 



< 



2« 111 0_ ~ 



<9p 



2s-l 



+ 



0,£ 



_9 
<9p 



+ 



o,s 



<9p 



u 



+ 

2s-2 



0,E 



hence by (I20I) . (|34|) . a trace theorem and the induction hypothesis, we have 



V, 



2s lll,n_ ~ 



< K 



2s-l, 



(2s - I)*" 1 + C u - 7 f_(2 S )! + C /+7 ^(2s)! < C(2 S )! 7 



2s 



for suitable constants C, 7 > independent of s. Therefore, from [3] we obtain for k > 2 



— K~ 



(73) 



(k-2 

U=o*" 



<9p 



2s-l 



+ 



<9p 



2> 



-» 



+ 

2s-2 



2s lli.n_ 



< 



1: 



2s-l 



+ " 



2s||^+2,0_ 



2s-2||£+2,n 4 



(2 S )! 7 



2s 
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Now, U2 S _ 2 = A( 2s )/ + (see eq. ( 120]) ). hence using (J7J) we have 



\u 



"2s-2\\i + i^ + 



= ll^ (2s, / + || w . n+ <E^ £ ll^ ,2 * , /■ 

4=1 ' \a\<i+l 
k k 

< E^E 7? +2s (l«l+ 2s)! < ^7^(^ + 23)! 

(2s)!(l + 7/ ) fc+2S 7r 



=1 \a\<l+l 



k+2s 



2s 



< fc 2 (2 S )£>5 

< (2s)!7 2s , 

for suitable 7 > independent of e. Also, by (152"]) we get 



< 



|^2s-l ( M 2s-1 M 2.s-l) l^il 



^ ||^2 S -l|| m) n_ + 



'2s-ilu+i,n. 



+ M 



2s-l|U+lo_ • 



Equation (|2"0|) gives = 0, and by (I2"0l) . (IMj) and the induction hypothesis, we obtain 



< 



2s-l 2s-l 



hi S 



(£ + l)!C ,m (2s-l) 2s ' 1 7 



Thus, (173!) becomes 



— Ivrl 



< ((£ + l)\C l+1 (2s - l) 2 -^ 2 - 1 + Cfi 1 + l)!(2s)! 7 2 i + (2s 



+C fe+1 (2s)! 7 2s 
< Cf +1 (2s) 2s 7 2s , 

for suitable constants Ci, 7l > independent of e. This establishes ( 168]) ; to establish (1691)- 
(170]) we will simply check that the boundary data in (162]) satisfies the appropriate bound (so 
that we may apply Lemma |3]). Since 

we see that for 6 G [0,L E ], (cf. fl62])) 



v£ (0, B) 



< \V 2 ~ (0,e)| + \u+ (O,0)| + |u*(0, 
C(2s) 2s 7 2s , 



< 

no 
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which is the bound that allows us to apply Lemma [3] and conclude that for V^a (p,0), the 
estimates ( 17D|) hold as desired. 

The case of even j: If j is even, then j : + 1 is odd and we would like to establish bounds 
for V 2 ~ +1 and V^s+i (with 2s + 1 = j + 1), which satisfy (1591) . ( 1601) respectively. First, for V 2 ~ +1 
we see from the variational formulation of fl59l) that 



|^2s+l||i if 



< 



< 



0,S 



C + (2sY s Y s < C + (2s + l) is+1 7 



2s+l,_,2s+l 



and, in a similar fashion as above, we obtain 



I k 2 1 1 1 ~ 

jfei l^2 S +i| fci n_ - ^ 1 X^lr 2 * 



£=0 
' fe-2 



f+l/2,£ 



fell 



< ^+1 J ^ (i±lH^+l (2s) 2, 7 2 S + c+(2s + 1)2 ,+1 7 2,4 



£=0 



,2s+l 



< ^ +1 (2 S + l) 2s+1 7r 

for suitable constants C 2 ,7 2 > independent of e. Finally, from the above result we see that 
the boundary data of (ISO}) satisfies the appropriate bound, hence by Lemma [31 V^ +1 satisfies 
f l69|) and (J2U as desired. ■ 

In view of the previous theorem, we define the (truncated) interface layer expansion(s) as 



(74) 
where 



M /L :— {yt,Mi v I,m) 



2M+1 



2M+1 



(75) ^ M = £ eiV+ , v' hM = ^ 

3=0 j=0 

The following corollary follows from Theorem HJ 



Corollary 5 There exist constants C, 7, G,X > depending only on the data, such that 
under the assumption e(2M + l)max{7,i^} < 1, the functions vj ~ M ,v} ~ M defined by (75\ ) 
satisfy 

\dldlvt^9)\<e^{p+lf/ 2 q \ f|Y, 
for p,q G N , p G [0, p s ] and G 5(G) given by 
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Proof. By (I75|) and Theorem [4] we have 

2M+1 2M+1 2M+1 

j=0 j=0 j=0 



and 



2M+1 2M+1 , v g 

i=o j=o ^ " ' 

/ 2 \ 9 2M+1 

< e-V(p + l) 1/2 g! J ^ (^( 2M + 



j=0 

2 N? 



Finally in this section, we wish to see what the contribution of the interface layers is, to 
the remainder of the expansion. For the interface layers in f2_ we easily see that 

(76) - A (vJ M ) + (vJ >M ) = 0. 

Now, by construction of the functions V£ we have (with the aid of (13 7p and (1431) ) 

oo 2A/+1 

Le Km) = E £i E 

i=2M+2 j=0 

2A/+1 oo 2M+1 oo 

= - E E : "i> 1 ' 1 Vy -EE £l p ~ 2 - J *i 2 - 3 d 2 e v; - 

j=0 i=2M+2 j=0 i=2M+3 

2M+1 oo 

-EE e^^deV;. 

j=0 i=2M+4 

By Lemma 2.12 of [8], we have the bound 

(77) \L e vt M (p, 6»)| < J^ M + 2 (e(2M + 2) + |p|) 2M+2 V (p, 61) G B Po (0) x 5(6), 

for some 6 > independent of s. (As before, B$(z) denotes the open disc in the complex 
plane of radius 5 centered at z, and S(Q) is given by (|4"91) .) 



Remark 2 Corollary [3] shows that the interface layer functions in Q + behave just like the 
boundary layers, while the interface layers in f2_ are smooth. This will be taken into consid- 
eration in the design of the approximation scheme in Section \4\ ahead. 
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3.4 Remainder estimates 



We now consider the remainder r £ = (rfj_,ri) in the decomposition ( !T0|) . which is given by 



(7? 



(79) 



r- = u - ~ w - ~ Xil v i,mi 



and by construction, satisfies the equivalent (but homogeneous) boundary conditions as u £ 
on dfl. To see this note that on <9f2_\£ we have 

( r -) = ( u - ~ w -- XilvJ, m ) | an _\ s = 0, 

by (@D, (E2D and (ESJ. On <9fi+\£ we have 

( r +) L + \£ = « - < - " Xjl^a/) | an+ \ E = 0, 

by ©, (J36D and flT^. Finally on £ we have 

« - r E _) | s = (t£ - ui - wi + + Xil^Im ~ Xil^i,m) | s = 
by ©, dSU) and 



<9z/ dv 



r.du 6 , du £ dw £ _ Q <9u>I dv r 

r — — h e — — H — 

<9z/ <9z/ <9z/ <9i/ <9f 



M 2 -f,M 
- £ - 



<9z/ 



by ©, and ( 153|) . Moreover, we have the following. 



Theorem 6 Let r e = (r+,rt) be given by $78\)-( 79) and let L £ r e + = — e 2 Ar £ + + r £ + and 



Lir e _ = —Art + r £ _. Then there exist constants K 1 ,K 2 > independent of e, such that 

(gQ) II T „e II <-/„/oi,f , o\ \2M+2 

and 



\L e r%\\ <(e(2M + 2)K 1 ] 



51) 



\Lir e _\\ on _<(e(2M + 2)K 2 



,2A/+2 



Proof. We first consider flSUl) and we have 



(82) 



L e r; = L £ (u% -w £ + - XblUbl ~ Xil^m) 

= L e (u% - w £ + ) - L £ (xbl u %l) ~ L e (Xil^im) 
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From ( 128]) we notice that 

(83) L £ {u%-wl)=e 2M+2 A^f +1 
and also 

L e (XblU £ bl ) = £ 2 ( a Xbl) u BL ~ 2z 2 Vxbl ■ Vw^ + Xbl L sU £ bl , 
where the function Xbl equals 1 for < p < p Q and for p > [p x + p )/2. Hence by (1471) . 

\\e 2 (A Xbl ) u% L \\ 0>n+ < e 2 (l + (e(2M + l)K) 2M+l ) e"^ 



and 



\s 2 Vxbl ■ V^kllo,^ + ^(2M + mf M+l ) e- ap/£ . 



,\ 2M+2 



for some appropriate constant K > 0. Therefore, by ( )48l) and the previous two inequalities, 
we obtain 

(84) ||L e (x B LU% L )\\ 0>n+ < (e(2M + 2)K) 2M+2 . 

In a completely analogous way, we may obtain bounds for L e (xil^m) > v ^ z - 

(85) (x/l^m) || 0in+ < (e(2M + 2)lf J 

for some appropriate constant K > 0. Combining (l82l) - (|85l) we have 

£ lk 2A/+2A(M+1) / + || ^ + + (e(2M + 2)A) 2M+2 + ( e (2M + 2)£) 2M+2 

< e 2M+2 7 Jf +2 (2M + 2)! + (e(2M + 2)A') 2M+2 + (e(2M + 2)kJ 

< (e(2M + 2)K l ) 2M+2 , 

for a suitable A'i > independent of e. This establishes (180j) . 
Turning our attention to (IHTj) . we have 

Lirl = Li (mI - wi - Xil v i,m) 

(86) = 
We have from (J55j) (with the aid of (J2SD) 

hence by (1761) . 

ll L i ri IU_ = ||Mx«^>)L_ 

= \\{&Xil) vJ tM - 2Vxil ■ Vf 7 " M + XiLLivJ tM \\ on 

< ||(Ax/L)^|| 0>n _ + \\2Vxil ■ VvJ jM \\ oa _ . 
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Since the function \il equals 1 for < p < p E and for p > (p 2 + p )/2 (cf. ([15]) ). we further 
get (using 



2M+1 2M+1 



|Lir £ || < eP (\\V~\\ + IIW~|| ^ < V 

i=o 

2M+1 

< Yl (ei(2M + l)) j <(eK 2 (2M + 2)) 



£ J j J ^f J 



x2M+2 
3=0 

for a suitable K 2 > independent of e. Thus ( IHTj) is established and this completes the proof. 



Remark 3 Theorem shows that for eM sufficiently small, the remainder in / TiOj) is expo- 
nentially small, hence it need not be approximated. This information will be utilized in the 
next section when we will construct the approximation to u e . 



4 Approximation results 

We begin this section with the variational formulation of ((TJ) — (JSJ) , which reads: Find u £ 
(u e + , u £ _) G Hq {9) such that 

(87) B £ (u e , v) = F(v) V v = «, v e _) G H] (Q) , 

where 

B £ (u £ ,v)= [ {e 2 V«)-V«)+«}+ f {V{u £ _)-V(v e _)+u £ _vt}, 



(89) F(v) = / f+v% + / f-vl_ + / hv. 

Jn+ Jn- Jt, 

It is straight forward to show that the bilinear form (188j) is coercive and continuous on Hq (Q), 
hence the variational problem (1871) admits a unique solution thanks to the Lax-Milgram 
lemma. The discrete version of (IHTj) reads: Find u £ N = (u^, u^[) G Vn C Hq (f2) such that 

(90) B £ (u%, v) = F(v) Vv = «, ul) G V N C ff a (fi) , 
and by Cea's Lemma we have 

(91) K-«Srll e < \\u £ -v\\ £ , 

v£V N 
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where the energy norm ||-|| is defined as 

(92) \\u\\ 2 E = B £ (u,u). 

We now describe the subspace Vn- For simplicity, we will focus on quadrilateral elements, 
even though triangular elements are also possible (see [7J for this and other choices of a suit- 
able mesh). Since the behavior of the solution u e depends on the value of e (cf. Theorem [1]), 
we distinguish between the cases upe > 1/2 and npe < 1/2 (with k61 a fixed constant) as 
follows: If npe > 1/2 then the mesh does not need any special design, as in this case the poly- 
nomial degree p of the approximating functions is high enough to ensure good approximability. 
Hence, in this case the mesh A only needs to be regular in the sence of j2] (or satisfy conditions 
M1-M3 in [7]). In the case npe < 1/2 the mesh will include elements of size 0(pe) along dQ + 
in order for the boundary and interface layer effects to be captured - these are referred to 
as needle elements in [7j. We now describe one such possible construction: Let be given 
by fll3p . and divide <9f2 + \X into subintervals (9j,9j+i) ,j = l, ...,m — 1,9 G dQ + . Then draw 
the inward normal at Oj of length p (see eq. (ITT]) ) and connect each point (p ■, 9j) = (p , 9j) 
using the curve p = p (=constant). Further, divide each 

(93) ) . := {(p, 9) : < p < p , 9 j < 9 < 9 j+1 } , j = 1, .., m 
into , (^+ 2 ) i , where 

(nf). = {(p,9):9 j <9<9 j+h 0<p<±p K P ey 

In the above definitions, k G 1 is a fixed constant, p is the degree of the approximating 
polynomials and we recall that we assume npe < 1/2. This will define a mesh 

A-= {(0^,(^1™, 
over Q°i_. We may define a completely analogous mesh over (see eq. (fT5|) ). as 

A|:={^') t .(^)JL. 

with (^s 1 )^ > (^s' 2 ) fc defined in an analogous way as Next, let {(nyj* 

be some subdivision of that is compatible with A^_ and A^, and define the mesh 

(94) a + = { (a% 1 ) k , (fi° E 2 ) fc , (n ; 1 )^ , (n°+ 2 ) . , , k = 1, n, j = 1, m , < = 1, e) , 

over Q + . The mesh A_ over f2_ is simply be chosen to be compatible with A + , and regular, 
in the sense of [2]. The mesh over the entire domain Q is then taken to be 

(95) A = A+UA_, 
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and we assume that the number of elements in A is bounded independently of s. The above 
mesh satisfies the definition of a regular admissible boundary layer mesh (Definition 3.2 in [7j), 
which implies the following: With S : = [0, 1] x [0, 1] the usual reference square, we associate 
with each quadrilateral ftf G A a different iable, bijective element mapping 



which, in this case, satisfies V i 



\\D a Mt\\ Loois) <i M \a\\\/aeN 2 . 
The space Vn is then defined as 

(96) V N = {u G H 1 (ft) : u| n ± = <f) p o (M±) _1 for p G Q P (S)} PI Hq (ft) , 

where Q P {S) denotes the space of all polynomials of degree p in each variable defined on the 
reference square S. Note that 

N = dim ^ = O (p 2 ) . 

Now, for p > 1 we define on the space of continuous function C ([0, 1]) , the operator ir p 
by interpolation in the p + 1 Gauss-Lobatto points, and on S we introduce the interpolation 
operator n p as the tensor product of the two one-dimensional operators n p and ir^. Then, by 
Lemma 3.8 in [TJ, we have that for any u G C°° (S) with ||D q m|| 05 < 7'°' |a|! V a G Nq, there 
exists a constant a > depending only on 7, such that 

(97) \\u - Tl p u\\ LOO{s) + || V (« - U p u)\\ LOC{s) < e-". 
Moreover, there holds (see, e.g., Lemma 3.7 in [7]), 



(98) 



H n p M llL-(S) ~ (1 + hip) 2 |M| LOO(5) 
|<9ji pw|| l°°(s) ' H^/Hp-uH L oo( S ) < p 2 (l + \np) 2 \\u\\ L oo(5) 



We now prove our main approximation result. 

Theorem 7 Let u £ G Hq (ft) , G Vjy 6e i/ie solutions of < f%5)) and (fP5|j, respectively, with 
Vn defined by If96]) on the mesh A <jwen 6?/ (TJ2J). Further, assume that dQ is analytic and the 
functions f± are analytic on ft± tu/w/e #ie function h is analytic on S. Then, for k sufficiently 
small, we have 

\\u £ -u N \\ £ <N 2 e^, 
for some constant b > independent of e and p. 
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Proof. We consider the cases npe > 1/2 (asymptotic case) and upe < 1/2 (pre-asymptotic 
case) separately. 

Case 1: Kpe > 1/2 (asymptotic case) 

By Theorem [1] there exist constants C, K > depending only on the data such that 

e\\D a u £ + \\ on+ + \\D a u £ _\\ Qn < CeK H max{\a\ ,£ _1 } H Va = 1,2 

Now, by Lemma 3.10 of [7J we have that for each element map Mf 



\D a KoM±)|| 0S < 7 L Q| |a|!e^, 



for some constants 7 ± > independent of e and i. Hence, by (J97l) and Lemma 3.6 of [7j we 
have 



|| (u% o ) - n p {u £ ± o Mf ) \\ Loo{s) + || V (« o Mf ) - n p («« o Mf )) || Loo(5) < e^ +1 / £ . 
Since 2/«p > we have —up + 1/e < — crp + 2/«p and, under the assumption that k < a/2, 



M e ± o Mf) - n p « o Mf) || ioo(g) + || V (« o Mf ) - n p {u% o Mf )) || Loc(5) < e"^, 

for some constant 6 > 0. By fl9Tl) and fl92l) we get the desired result. 
Case 2: Kpe < 1/2 (pre-asymptotic case) 

In this case we utilize the expansion and regularity results of Section [3] which state that 
the solution u £ = (u+, uf_) can be written as 

U % = W M + XBLUBL + Xl L vi,M +r +> 
U £ _ = W M + Xl L v7,M + r -i 

with each term defined and analyzed in subsections I3.1H3.41 We begin be selecting M in such 
a way that eM is sufficiently small for all the regularity results of subsections 13 . 1HTT41 to hold 
true. (The lack of concreteness on our part is due to the careful constant selection made in 
[7], i.e. such a choice for M is possible by jJJ). The proof relies on the following observations: 

1. The terms w^, vJ M are analytic in their respective domains, hence (197|) may be applied. 

2. The construction of the mesh allows us to approximate u^ L and vf M at an exponential 
rate. 

3. The choice of M renders the term ||r e || e ^ negligible (exponentially small), hence the 
remainder need not be approximated. 
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Let us first consider item 1 above. For the term v T M we have, by Corollary [H that for 
each element map Mr 



\D a fc/°MT)|| 05 <C H 



all, 



hence by ( T97f) . 



vl M ° Mr) - n p {vl M o M~) \\ loo(s) + || V {{vy M o Mr) - U p (vJ M o M~)) || ioo(s) < e"*. 



The same works for the other two terms (the details are ommitted). 

Next let us comment on item 2; since the steps are very similar for both Ug L and vf M , we 
will only consider the latter. Without loss of generality we assume that Xil * s 1 i n an d 
otherwise. Hence we only need to approximate vf M within f^. To this end, we note that 
the mesh in consists of two types of elements (cf. ((94 



{Qf). = |(p, 9) : 6, < 6 < J+1 , 0<p< \ Po ^ and (flj 2 ) . = ) . \ (fi /) . . 

For (fi^' 1 ) • ' w ith associated mapping M ' , we have from Proposition 3.11 in [7] 

\\D a Km o r X o Mf) || ioo(5) < e^l Q l |a|! Va G Ng, 
where ^ was defined by (|T2|) . Therefore, by (1971) 

II o v;- 1 o m;* 1 ) - n p o ^T 1 o m;- 1 ) || loc(s) + 

+ || V ((vf M o r 1 o Mf) - n p {vf jM o V;" 1 o Mf)) \\ Lao(s) < e "* e -* 

from which the desired result follows provided n < b. Now, let us consider the approximation 
of vf M over the elements (f^ 2 ) , with associated mapping Mr . From ( )47|) we have 

IIx/l^m ° ^ 1 o M^\\ Loo{s) < C a e~^ , || V (xiiftjt ° iT 1 ° M D IL (5) £ O^-V'"*. 
Therefore, from (j98l) we get 

IIXil^m ° ^ o M;' 1 || ioo(5) < (1 + lnp) 2 e-^ 

and 

II V {xilKm ° ° M^ 1 ) || ioo(5) < e-y (1 + hip) 2 e "* 
from which the desired result follows once we use ( 1921) . 

Turning to item 3, we have from the variational formulation (IH71)-(l89l). that the remainder 
r £ = {r e +) r e _) satisfies 

B £ (r £ ,v) = / L £ r £ + vdx+ / L x r £ _vdx, Vi> G Hq(Q). 
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Hence by Cauchy-Schwarz's inequality we get 

S £ (r E ,r £ )<||L £ r;||^ + + 11^111^, 
and by Theorem [6] we obtain 

||r £ || £ = [B e (r £ ,r £ )} 1/2 < (e(2M + 2)K) 2M+2 < e~ bp , 
for some suitable constant b > 0, depending only on the data. This completes the proof. ■ 

5 Numerical results 

In this section we will illustrate our theoretical findings for the model problem dl])-©, in the 
case when / + = /_ = l,h — and the domain Q consists of the two subdomains Q + and 
delimited by the three concentric circles with radii 1, 2 and 3. In other words, Q + is the 
domain inside the two concentric circles of radii 1 and 2, while f2_ is the domain inside the 
two concentric circles of radii 2 and 3, as shown in figure 2. 




Figure 2: Domains f2+ and f2_ used for the computations. 

We expect to have a boundary layer along <9f2 + \X (the circle of radius 1) and an interface 
layer along £ (the circle of radius 2). The mesh, shown in figure 3, accounts for the presence of 
the layers by including thin elements of size pe along dQ + \E and £ - the value of the constant 
k appearing in the definition of the mesh in the previous section was taken to be 1 (a value 
known to produce almost the same results as those obtained with the "optimal " value of k, 
see, e.g., [H]). An exact solution is available for this problem, hence our computations are 
reliable. 

The computations were performed with the commercial package StressCheck (E.S.R.D., 
St. Louis, MO) which is a p-version FEM software package allowing the polynomial degree to 
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Figure 3: Design of the mesh (on a quarter of the domain). 




Figure 4: Approximate solution for p = 8, e = 0.01. 

vary from p = 1 to p = 8 (on a fixed mesh). Figure 4 below shows the approximate solution 
for p — 8, £ = 0.01, and figure 5 shows the convergence (in the energy norm) as p is increased 
- the exponential convergence is readily visible. 

6 Conclusions 

We have studied the finite element approximation of a singularly perturbed transmission 
problem posed on a (smooth) domain with analytic boundary. Upon obtaining appropriate 
regularity results, via asymptotic expansions, we were able to design and analyze an hp finite 
element method for the robust approximation of the solution to the singularly perturbed 
transmission problem. We showed that under the assumption of analytic data, our method 
converges at an exponential rate, independently of the singular perturbation parameter. This 
is in line with our one- dimensional results [12], as well as with two-dimensional results for non- 
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Convergence fore = 10 ' 

10 2 , , , i , i i ,— | , , i i , i I,, Convergence for e = 10" 




Figure 5: Convergence of the approximate solution: loglog plot (left); semilog plot (right). 



transmission problems [7] . The approximation of singularly perturbed transmission problems 
on non-smooth domains is the focus of our current research efforts. 
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